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The Problem 


How do we generate a complex time-evolving sound composed 
of many discrete components (ie. harmonic, or inharmonic, 
partials)? 


The brute force approach: 

Use one sinewave oscillator plus a pair of envelopes 
(amp, freq) per partial, then mix all the sources 
together. 

The elegant solution: 

Find a way of combining a few simple sources (ie. 
sinewave oscillators) to generate lots of components 
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The Problem, mathematically 

stated 


N -1 

E a*cos(u) k +<|> k ) 

k= 0 

In other words, we want to generate a sound with N components 
summed (mixed) together (2), each component with: 
amplitude a k 

frequency co k 
phase (j) k 


The simplest case is the Fourier series of a pulse 
with N harmonic components 


N 

Z cos (k U3 0 ) 

k =1 
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Closed-form Summation 


The first of these ways to combine sinusoids to create a harmonic 
series is to employ a well known mathematical device: a closed- 
form sum. This is a well known one: 




1 —r 


2N + 1 


1— r 


For the harmonic series, we have a well-defined expression: 
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Discrete Summation Formulae 


Moorer provided more general principles for closed-form 
sums. In particular, he came up with several expressions that 
could generate a variety of spectral combinations, for instance: 


sin (to) — g sin (co — 0) 
1—2gcos(0)+g 2 


^ g k sin(co + k0) 

0 


In general, for a bandlimited spectrum, we have 


sin((jo) —gsin((jo —0) —^ iV sin((jo + [iV + l]0)+g iV+1 sin((jo + iV 0) 

l-2gcos(0) + g 2 


N-1 

^ g k sin(uo-hk6 ) 

k = 0 
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Implementation Example 

opcode NBlsum,a,kkkki 

ka,kw,kt,kk,itb xin 

aphw phasor kw 

apht phasor kt 

asinl tablei aphw,itb, 1 

asin2 tablei aphw - apht,itb, 1,0,1 

acos tablei apht,itb, 1,0.25,1 

ksq = kk*kk 

asig = (asinl - kk*asin2)/(l - 2*kk*acos + ksq) 
knorm = sqrt(l-ksq) 
xout asig*ka*knorm 

endop 
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Waveshaping 

Waveshaping, or non-linear mapping, is a method of producing 
lots of components by distorting the shape of a sinusoid with a 
function: 

f( cos(co)) 

if f(x) is non-linear (ie. not a straight line), then a number of 
harmonic components will be produced. For instance: 



l -gx 
l-gx + g 2 


f( cos(to)) 


1 — gcos(io) 

1 —2gcos(co)+g 2 


= X 9 k cos (k w) 

k =0 
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Hyperbolic Tangent 
Waveshaping 

Another distortion method that can be useful for creating classic 
analogue waveforms is the hyperbolic waveshaping method: 


square (w)=tanh(ksin(co m )) 


By limiting the amount of modulation, it is possible to produce 
nearly bandlimited square waves. From these it is simple to 
produce sawtooths: 

saw (<x>)=square{ to)(cos (to)+l) 
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Implementation Example 


opcode Waveshape,a,kkkiii 
kamp,kf,kndx,isin,itf,igf xin 
asm oscili 0.5*kndx,kf,isin 
awsh tablei asin,itf, 1,0.5 
kscl tablei kndx,igf,l 
xout awsh*kamp*kscl 
endop 



; function tables: 
f2 0 16385 "tanh” -157 
157 

f3 0 8193 4 2 1 


opcode Sawtooth,a,kkkiii 
kamp,kf,kndx,isin,itf,igf xin 
amod oscili l,kf,l, 0.25 
asq Waveshape kamp*0.5,kf,kndx,isin,itf,igf 
xout asq*(amod + 1) 
endop 
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Asymmetrical FM 

Palamin et al built on previous work by Moorer to propose this 
interesting formula for changing the symmetry of the double¬ 
sided FM spectrum 


1 1 

s(t) = exp(0.5k(r—) cos (co m )) sin (co c +0.5k (r + —)sin(co m )) 

r r 


oo 


= Z r n j„(k)sin(u) c +7iu3 m ) 


n =—co 


As the expanded sum demonstrates, this includes a new 
scaling variable r, that will shift the symmetry of the 
spectrum away from the carrier frequency. 
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Implementation Example 

opcode Asfm,a,kkkkkii 

kamp,kfc,kfm,knx,kR,ifn,ifn2 xin 

kndx = knx*(kR+l/kR)*0.5 
kndx2 = knx*(kR-l/kR)*0.5 

afm oscili kndx/(2*$M_PI),kfm,ifn 
aph phasor kfc 
afc tablei aph+afm,ifn, 1,0,1 


amod oscili kndx2, kfm, ifn, 0.25 
aexp tablei (abs(kndx2) - amod)/50, ifn2, 1 
xout kamp*afc*aexp 

endop 


The exponential is implemented with a function table 
lookup. In order to use the table lookup limiting mechanism, 
we draw up an inverse exponential from 0 to -50 and then 
reverse the sign of the argument to it. 


;(inverse) exponential function table 
f5 0 131072 "exp" 0 -50 1 
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Phase Aligned Formant Synthesis 


Puckette, in 95, proposed PAF as ring-modulation of a sinusoid 
with a complex wave. This is actually similar to one of 
Moorer’s double-sided DSF equations 


00 


00 


^ g |/c| cos((jo c +k(jo 0 )=cos((jo c )[l + l^ g k cos(k<x> 0 )] 


k =—go 


= cos (co )[ 


±-9 


l-2gcos((« 0 ) + g‘ 


k = 1 

] = c(iu )M(io 0 ) 


The complex modulating wave can be implemented using 
waveshaping, so this is a distortion technique. For this 
purpose the formula for M(.) above is rearranged into: 


M(oo) = r±£f(2 % /g Sm ^ /2 ^ ), with f(x) = 


i -g 


i -g 


l + x‘ 
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Implementation Example 


opcode PAF,a,kkkkki 


kamp,kfo,kfc,kfsh,kbw,itb xin 

kn = int(kfc/kfo) 

ka = (kfc - kfsh - kn*kfo)/kfo 

kg = exp(-kfo/kbw) 

afsh phasor kfsh 

aphs phasor kfo/2 


opcode Func,a,a 


asig xin 



al tablei 2*aphs*kn+afsh, 1,1,0.25,1 
a2 tablei 2*aphs*(kn+l)+afsh, 1,1,0.25,1 
asin tablei aphs, 1, 1, 0, 1 
amod Func 2*sqrt(kg)*asin/(l-kg) 
kscl = (l+kg)/(l-kg) 
acar = ka*a2+(l-ka)*al 
asig = kscl*amod*acar 
xout asig*kamp 

endop 
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Modified FM Synthesis 

The modified FM synthesis method is a variation on classic 
FM, which exhibits modified Bessel Functions in its expansion: 

oo 

e kcos< “")- k cos(co c ) = e- k ^ / n (k)cos(io c +nco m ) 


Its advantage is that modified Bessels are unipolar and 
decreasing, which allows for a better-behaved spectrum, with 
similar computational costs to FM (ie. very little). 

Applications for bandlimited classic analogue waveform 
generation, among other uses, have been proposed. 
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Implementation Example 


opcode ModFM,a,kkkkii 


kamp,kfc,kfm,kndx,isin,iexp xin 
acar oscili kamp,kfc,isin,0.25 
acos oscili l,kfm,isin,0.25 
amod table -kndx*(acos-l)/50,iexp,l 
xout acar*amod 



;(inverse) exponential function 
table 

f5 0 131072 "exp" 0 -50 1 


endop 
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Conclusion 


Distortion techniques provide efficient and elegant ways of 
synthesising sounds. With the advent of other methods, 
research interest in these had vanished for a number of years. 
Some novel formulations and applications have rekindled 
interest in the area. 

Although the basic techniques flourished in the brave days of 
early computer music, there are a number of interesting 
possibilities still left to be pursued. 
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